Accepted by the Astrophysical Journal 

Preprint typeset using I^T^X style emulateapj v. 5/2/1 1 



SLOSHING COLD FRONTS IN GALAXY GROUPS AND THEIR PERTURBING DISK GALAXIES: AN X-RAY, 

OPTICAL AND RADIO CASE STUDY 

Fabio Gastaldello^'^, Laura Di Gesu^ '^ Simona Ghizzardi^ Simona Giacintucci^^ Marisa Girardi^'^, Elke 

ROEDIGER^'^, MARIACHIARA ROSSETTI^'^ FABRIZIO BRIGHENTI^^'^\ DAVID A. BUOTE^, DOMINIQUE ECKERT^^, STEFANO 

Ettori^^'^^, Philip J. Humphrey^, and William G. Mathews^^ 

Accepted by the Astrophysical Journal 



en 

o 

<N 
< 



o 
u 

o 

<N 
> 

00 

in 

o 
m 



% 



ABSTRACT 

We present a combined X-ray, optical, and radio analysis of the galaxy group IC 1 860 using the currently 
available Chandra and XMM data, literature multi-object spectroscopy data and GM7?r data. The Chandra and 
XMM imaging and spectroscopy reveal two surface brightness discontinuities at 45 and 76 kpc shown to be 
consistent with a pair of cold fronts. These features are interpreted as due to sloshing of the central gas induced 
by an off-axis minor merger with a perturber. This scenario is further supported by the presence of a peculiar 
velocity of the central galaxy IC 1 860 and the identification of a possible perturber in the optically disturbed 
spiral galaxy IC 1859. The identification of the perturber is consistent with the comparison with numerical 
simulations of sloshing. The GMRT observation at 325 MHz shows faint, extended radio emission contained 
within the inner cold front, as seen in some galaxy clusters hosting diffuse radio mini-halos. However, unlike 
mini-halos, no particle reacceleration is needed to explain the extended radio emission, which is consistent with 
aged radio plasma redistributed by the sloshing. There is strong analogy of the X-ray and optical phenomenol- 
ogy of the IC 1860 group with two other groups, NGC 5044 and NGC 5846, showing cold fronts. The evidence 
presented in this paper is among the strongest supporting the currently favored model of cold-front formation 
in relaxed objects and establishes the group scale as a chief environment to study this phenomenon. 

Subject headings: gSilsixies'. clusters: intracluster medium — galaxies: groups: individual (IC 1860, NGC 5044, 
NGC 5846) — hydrodynamics — X-rays: galaxies: clusters 



1. INTRODUCTION 

The current generation of X-ray telescopes (XMM and to 
a greater degree Chandra) with their angular resolution and 
sensitivity has revealed a wealth of small scale features in the 
intra-cluster (ICM) and intra-group medium (IGM). One par- 
ticularly interesting is the phenomenon of the "cold fronts", 
i. e. sharp surface brightness discontinuities, interpreted as 
contact edges between regions of gas with different entropies 
(see the review by Markevitch & Vikhlinin 2007]) . Cold fronts 
appear to be almost ubiquitous in galaxy clusters ( [Markevitch 



eraLl|2003t [Ghizzardi et~aLl[20!0l [Lagana et aL][2010| ) and 
they are found both in dynamically active objects and in cool 
core relaxed clusters. In merging clusters cold fronts arise 
during merger events through ram-pressure stripping mecha- 
nisms that induce the discontinuity among the merging dense 
subcluster a nd the less dense surrounding ICM (e.g. Vikhl inin| 
et al.|200T ). In relaxed clusters cold fronts are most likely m- 
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duced by minor mergers that produce a disturbance on the 
gas in the core of the main cluster, displace it from the cen- 
ter of the potential well, and decouple it from the underly- 
ing dark matter halo through ram pressu re (Markevitch e t aL| 
[200Tl|Churazov et al.||2003t|Ascasiba7& Markevitch..2005r 
The oscillation of the gas of the core around the minimum 
of the potential generates a succession of radially propagat- 
ing cold fronts, appearing as concentric edges in the surface 
brightness distribution of the cluster. These fronts may form 
a spiral structure when the sloshing direction is near the plane 
of the sky. When the sloshing direction is not in the plane of 
the sky concentric arcs are observed. The sequence of events 
is described in great detail in the simulations presented in |As-| 
Icasibar & Markevitch] ( |2QQ6D . 

Cold fronts are unique t ools to understand the physica l 
properties of the ICM (e.g. Markevitch & Vikhlinin 2007|. 
They ca n also be in principle used as a gauge of the merger 
activity ( [Qwers et al.||2QQ9a| ), in particular of the more fre- 
quent minor merging of smaller subsystems compared to the 
less frequent and more spectacular major mergers. Cold fronts 
do not need a favorable geometry for being detected as the 
more elusive shock fronts (e.g. Markevitch 2010) and they 
may provide the only indication at X-ray wavelengths of a 
minor merger in a seemingly relaxed cluster. Cold fronts can 
be used, with the help of simulations, to infer many merger 
characteristics, as for example from the presence and size of a 
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spiral pattern in the surface brightness distribution t he direc- 
tion and time of peri-center passage of the perturber ([Ascasi- 



2QTTr 



bar & Markevitch|2006] [Johnson et al.|201Q|[Roediger et al. 



Recently, gas sloshing in the core has also been invoked to 
explain the formation of diffuse radio minihalos in relaxed, 
cool-core clusters. These faint, steep- spectrum radio sources 
are relati vely rare, with only few clusters with confirmed de- 
tections ( [Feretti et al.||2Q12| . A spatial correlation between 
the minihalo emiss ion and cold fronts has bee n observed in a 
few clusters ( Mazzotta & Giacintucci 2008 ; Gi acintucci et al.| 
[2013), with the minihalos contained within the region con- 
nned by the cold fronts, suggesting a tight connection be- 
tween sloshing motio ns and the origin of minihalos. Recently, 
ZuHone et al. p013 ) showed with high resolution MHD sim- 
ulations that gas sloshing can lead to turbulent reacceleration 
of relativistic electron seeds (e.g., from past AGN activity) 
and produce radio diffuse radio emission within the envelope 
of the sloshing cold fronts. 

Another complementary path to the study of presence and 
dynamics of cluster mergers uses optical data to perform a 
spatial and kinematical analysis of member galaxies. These 
types of studies allow us to reveal and measure the amount of 
substructure, and to detect and an alyze possible pre-m erging 
clumps or merger remnants (e.g., [Girardi & Biviano|[2002, ). 
This is certainly interesting for major merger clusters which 
are very rich in optical substructures (see for example the 
DARC - Dynamical Analysis of Radio Clusters - project, [GT] 
|rardi et al. 2007 ) but even more in the case of the minor merger 
scenario for sloshing. The assumption of a single Gaussian 
distribution is a good description for the galaxy velocities of 
the main cluster and the recently merged perturber may be 
ident ified by means of o ptical/dynamical substructure search 
(e.g., |Qwers et al.|2 011a). Optical studies can also investigate 
the presence of peculiar velocities of the bright central galaxy 
(BCG). In fact one o f the observable effects of the m erger 
scenario discussed by [Ascasibar & Markevitch| ( |2006| ) is the 
fact that, if the BCG sits at the peak of the dark matter dis- 
tribution, it is expected to start oscillating along with it after 
each subcluster flyby. Gas sloshing and BCG peculiar veloc- 
ities are caused by the same minor mergers. Given the in- 
creasing number of objects showing sloshing cold fronts this 
can explain also the large and apparently puzzling numb er of 
systems showi ng peculi ar velocities of thei r BCGs (e.g., |van| 
|den Bosch eT al. 2005i iCoziol et al.||2QQ9l ). As a matteroF 
fact Miralda-Escude] ( |1995| ) looking at the early evidence of 
BCGs spatially coincident with the centers of dark matter ha- 
los producing lensing features, commented that "the fact that 
cD galaxies often have large peculiar velocities relative to the 
average of the cluster galaxies has been used as an argument 
against their identification as cluster centers: however, clus- 
ters are continuously merging, and their density peaks do not 
need to coincide with their centers of mass. Substructure will 
cause the density peaks to move, in response to the gravita- 
tional forces of the in-falling material." 

It is yet another manifestation of the hierarchical nature 
of structure formation the fact that sloshing cold fronts have 
been detect ed also at the smaller mass scales of poor clusters 
(e.g., Virgo, Simion escu et al.|2010) and of groups of galaxies 
dGastaldelloetal. 12 009; Randall et al. 2009; Machacek eTaL] 
201 1| ). At this mass scale even a single massive galaxy can 



The aim of this paper is to investigate sloshing features, 
such as cold fronts and spirals in the surface brightness distri- 
bution, in the nearby. X-ray bright galaxy group IC 1 860. We 
will look for other observable properties of the sloshing sce- 
nario such as a peculiar velocity of the BCG and extended ra- 
dio emission possibly related to the presence of sloshing. We 
present in detail the results of the available XMM and Chan- 
dra data of the object, the available optical data in the litera- 
ture and new Giant Metrewave Radio Telescope (GMRT) data. 
We will then compare the observed features in IC 1860 with 
two other groups hosting cold fronts, NGC 5044 and NGC 
5846. 

The group IC 1860 has been originally classified as a 
Dressier cluster (DC-02 47-31, [Dre ssier 1980) and it is also 
known as Abell S301 ( Abell et al.j|1989| ). It appears then 
in many other catalogues, in particular the catalogue of |Maia 
et al. (1989). It is a nearby group of galaxies at z=0.022 (at 
a distance of 97 Mpc with the cosmology adopted in this pa- 
per) and its optical structure/membership has been studied in 
[Dressier & Shectman| ( |1988j ) and in particular, given the many 
additional redshift^ rovided by the 2dF redshift survey, by 
Burgett et al.|([2004|. It has been dis covered in the X-r ay band 



by the Einstein X-ray Observatory ( [Burstein et al.|1997) and 
ROSAT showed the presence of extended a nd diffuse group 
X-ray emission at a temperature of ^1 keV ( [Mulchaey et al. 
2003| ). Chandra imaging showed no evidence of AGN in 



duced activity in the form of cavities in the core ([Dunn et al. 



[2010| ). In our previous XMM analysis ( [Gastaldello et al.[2007| ) 
we mainly focused on the azimuthally averaged radial prop- 
erties of the system for the purpose of mass analysis, however 
noticing the presence of a possible disturbance in the surface 
brightness distribution. 

The paper is organized as follows: in ^we summarize the 
XMM and Chandra data reduction and analysis; in ^we an- 
alyze the X-ray images and surface brightness profiles finding 
evidence of a set of two surface brightness discontinuities and 
an excess brightness spiral feature; in ^we perform a spec- 
tral analysis of the temperature structure of the group finding 
evidence that the surface brightness discontinuities are cold 
fronts; in ^we model the pressure jumps between the cold 
fronts and m ^ we discuss the possible systematics of our 
X-ray analysis; m ^and in ^we present the analysis of the 
available optical ana GM7?r data, respectively. In ^we use 
simulations to constrain the merger geometry of IC 1 860 and 
in ^TQ| we discuss our results, making comparisons with the 
groups NGC 5044 and NGC 5846. Additional analysis per- 



formed on these groups is presented in appendix|A| in §11 



be the responsible perturber and the X-ray signature of slosh- 
ing cold front is a signature of the gravitational interaction 
between the galaxy (and its halo) and its host group. 



we 
summarize our results. We also report on the narrow~Sigle 
tail radio galaxy IC 1858 in appendix [C] 

The cosmology adopted in this paper assumes a flat uni- 
verse with Ho= 70 km s"^ Mpc"\ (^ni= 0.3 and (1a= 0.7. All 
the errors quoted are at the 68% confidence limit. At the dis- 
tance of IC 1860 V corresponds to 26.7 kpc. 

2. OBSERVATIONS AND DATA PREPARATION 

IC 1860 has been observed by XMM on February 04th 2003 
(obsID 0146510401) with the EPIC MOS and pn cameras for 
39 and 38 ks respectively and by Chandra with the ACIS-S 
configuration on September 2009 (obsID 10537) for 40 ks. 
We used both Chandra and XMM available archival data and 
we took advantage of the different characteristic of the two 
satellites. We exploited the wider field of view provided by 
XMM to analyze the source structure at larger radii while the 
higher resolution of Chandra allowed us to explore the de- 
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Figure 1. Exposure-corrected mosaic of the XMM MOS 1 and MOS2 images 
Gaussian smoothed on a 3" scale. Superposed over the image are the angular 
sectors used for the surface brightness profiles of Fig. [2]and discussed in the 
text. 

tails of the source in the inner region. \Ve processed the data 
following the standard SAsRand CIAOJ^hreads; in the fol- 
lowing subsections we briefly describe the data preparation. 

2.1. XMM 

Starting from the Observation Data Files retrieved from 
the XMM-Newton archive, we generated calibrated event files 
with SAS, version 1 1.0, using the tasks emchain and epchain. 
We removed proton flares both from the hard and the soft 
light-curve, using respectively a fixed th reshold criterium and 
a cr-clipping technique, as described in |Rossetti & Molendi 
([2010). We also checked for residual contamination by a qui- 
escent soft-proto ns component, ev aluating the in-over-out ra- 
tio as defined in 'De L uca & Molendi| ( |2004|) . Net exposure 
time after protons flare correction is 34.5 ks tor MOS-1, 34.8 
ks for MOS-2 and 29.4 ks for pn. We checked the observa- 
tion for contamination by solar wind charge exchange: the 
ACE (Advanced Composition Explorer) SWEPAlVp^proton 
flux was less than 4 x 10^ protons s~^ cm~^ and the SWICS 
O^^ 10^^ ratio was l ess than 0.3, values which are typical of 



the quiescent Sun ( Snowden et aL||2004| ). The spectra of the 
out-of-field-of-view events of CCD 4 of MOSl showed an 



anomalously high flux in the soft band (see |Kuntz & Snow- 
den|2008 ) and it was therefore excluded from our analysis. 

We then filtered the events according to standard pattern 
and flag criteria and we performed out-of-time correction for 
pn. For each detector we created images in the 0.5-2 keV 
band and exposure maps. We detected point sources using 
the task ewavelet and masked them using a circular region, 
centered at the source position and with a 25 '^ radius. Using 
the task emosaic, we combined the MOS images into a single 
exposure-corrected image shown Fig[T] and discussed in Sec 

2.2. Chandra 

In order to ensure the most up-to-date calibration, we repro- 
cessed all Chandra data starting from level 1 event files, us- 
ing the X-ray analysis package CIAO 4.3 in conjunction with 

^^ http://xmm.esa.int/sas/current/documentation/threads/ 
^ ^ http : //cxc .harvard . edu/ciao/threads/index . html 



the Chandra calibration database (Caldb) version 4.1.3. We 
took into account time-dependent drift in the detector gain and 
charge transfer inefficiency as implemented in the CIAO tools. 
To clean the data from periods of enhanced background, we 
performed a a-clipping to the light-curve in the 0.5-7 keV 
band. The final exposure time after deflaring is 35.7 ks. Point 
sources were detected with the the CIAO tool wavdetect and 
then excluded from all the following steps of the data anal- 
ysis using an appropriate mask region. We used the blank- 
sky dataset provided in the CALDB to perform background 
analysis, after properly re-processing and re-projecting it. To 
take into account possible temporal variations of instrumental 
background, we rescaled the background file for the the ratio 
between the count-rate of the observation and the background. 
To calculate the proper rescaling factor, we extracted spectra 
for both background and source files from an external region, 
not contaminated by source emission and we quantified the 
count-rate ratio in the hard band 9-12 keV. We then generated 
an image the 0.7-2.0 keV band and corrected it for its expo- 
sure map. The analys is we performed on the Chandra image 
is discussed in section [T2l 



3. X-RAY IMAGES AND SURFACE BRIGHTNESS 
PROFILES 

The sloshing scenario has been introduced to explain a 
number of features present in the surface brightness distri- 
bution of relaxed clusters and groups, in particular the pres- 
ence of concentric sharp surface brightness edges which have 
a temperature jump consistent with being cold fronts and the 
presence of an excess corresponding to spiral- or arc-shaped 
brighter regions inside the actual fronts. The morphology de- 
pends mainly on the angle between our line-of-sight (LOS) 
and the orbital plane of the subcluster: a spiral like structure 
is seen if the interaction is seen face-on, if on the contrary the 
line of sight is parallel to the orbital plane arcs are seen on 
alternating sides of th e cluster core ( Ascasibar & Markevitch 



^^ the data available at http://www.srl.caltech.edu/ACE/ACS/ 



2006[|Roediger et aLpOl 1 ). We will investigate in the follow- 
mg sections the surface brightness distribution of the X-ray 
emission of IC 1 860 as revealed by the XMM and Chandra 
data to look for such features. 

3.1. The outer surface brightness edge 

Visual inspection of the mosaic image of the two MOS cam- 
eras (Fig[lJ clearly reveals the presence of a sharp brightness 
discontinuity in the north-west direction (RA. 250° -340°, po- 
sition angles measured from the N direction). In the oppo- 
site direction (RA. 70° -160°), at large radii the group is more 
elongated and in the inner region a luminous plume is seen: 
this feature is sharper in the Chandra image so we will focus 
on it in section [32] To confirm this qualitative impression, we 
compared the surface brightness profiles in four angular sec- 
tor 90° wide as depicted in Fig[T] An evident slope change is 
seen in the 250°-340° profile, at a radius of - 170'' (Fig. [2]); 
conversely the profiles of the other three control sectors are 
fairly smooth with a continuous derivative. In the right panel 
of Fig. |2] we provide a zoom over the interesting radial region 
where tne edge is seen. The feature is also detected in the pn 
and the ACIS data (Fig [3]): however it appear less evident in 
the Chandra image, being close to the outer edge of the CCD. 

3.2. The inner surface brightness edge 

The Chandra image of Fig |3] provides a more detailed in- 
sight over the inner region oiIC 1860, confirming what al- 
ready pointed out in section |3.1[ in a narrow sector at R A. 
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Figure 2. Left panel: XMM surface brightness profiles in the selected angular sectors discussed in the text (section j3.1) and shown in Fig. [T] Green circles: 
sector P. A. 340° -70°. Blue circles: Sector P. A. 70° -160°. Black circles: sector P. A. 160° -250°. Red circles: sector!^. 250° -340°. The outer edge can be 
easily seen in the sector P. A. 250° -340°. Right panel: Zoom over the interesting radial region for the outer edge. The slope of the source profile is steeper in the 
sector PA. 250° -340° than in the control sector PA. 340° -70°. 

and seen in many hyd rodynamic simulation (e.g., the snapshot 
at 1.9 Gyr of Fig. 7 of Ascasibar & Markevi tch 2006) seen al- 
most face on. The southern edge of this "mushroom head" 
coincides with the inner surface brightness discontinuity. 



4. X-RAY SPECTRAL ANALYSIS 

To further investigate the nature of the brightness discon- 
tinuities we analyzed the temperature profiles across them. 
For this purpose we extracted XMM and Chandra spectra in 
several radial binq^ along the angular sectors containing the 
edges (RA. 250° -340° and RA. 130° -160°) and in other two 
control sector (RA. 340° -70°; RA. 160°-250°) in order to 
measure the temperature trend across the edges and check if it 
is actually different from the overall undisturbed temperature 
profile. 

We extracted XMM spectra generating a response file and an 
ancillary response file using the standard SAS tasks rmfgen 
and arfgen in extended source mode; we extracted Chandra 
spectra generating count- weighted spectral response matrices 
appropriate for each region using the task specextract. For 
each region, we extracted also background spectra from the 
blank-field event files. All spectra were re-binned to ensure a 
minimum 20 counts bin~^ and fitted in the 0.5-5.0 keV band. 




Figure 3. Chandra 0.7-2.0 keV image Gaussian smoothed on a ?>" scale 
and corrected for exposure map. Superposed over the image are the angular 
sectors used for the surface brightness profiles of Fig. |4]and discussed in the 
text. 

130° -160° the source is elongated in a luminous tail ending 
in a surface brightness edge. In Fig|4]we compare the profile 
of the source in this sector, with the profile of the source in 
other two control sectors as depicted in Fig|4] The profile of 
the source in the sector at RA.130°-160°is clearly different 
with respect to the other control profiles in the radial range 
20''-100''; the feature is more evident in the zoom of the left 
panel of Fig|4j a peak correspondent to the luminous tail is 
seen, then at a radius of ~ 100'^ the profile flattens out. The 
feature is clearly seen also in the XMM data. 

3.3. The spiral feature 

In Fig|5]we present the surface brightness residual map be- 
tween the XMM data and the best-fit radially symmetric two- 
dimensional /3-model describing the surface brightness distri- 
bution obtained with Sherpa: the residual map highlights the 
appearance of sloshing features. A characteristic spiral pat- 
tern in surface brightness can be seen, with the tail being con- 
fined by the outer cold front. The qualitative impression is that 
the overall morphology is tilted with respect to the plane of the 
sky, with the mushroom structure formed by the coolest gas 



with an APEC ( [Smith et al.||2001 ) plasma model as imple- 
mented in Xspec version 11.0. We take into account galactic 
absorption using a phahs model with the hydrogen column 
density frozen to the value provided in Kalberla et al.|( [2005| ). 
In section[6]we will explore the sensitivity of the results to the 
above assumptions. 

In Table [T] we present the results of the spectral fits in the 
interesting radial bins inside and outside the edges and in the 
same radial bins in the control sectors. In Fig[6]we show the 
temperature profile along the sectors containing the edges as 
obtained with XMM (left panel) and Chandra (right panel): 
the vertical dashed lines mark the radius of the two brightness 
discontinuity. The region outside the outer edge is covered 
only by the larger XMM field of view; therefore the complete 

^^ The range of the radial bins used for the spectral analysis are: 20^^ -40^^ 
40''-60'^ 60''-90'^ 100''-150'^ 180''-230'^ 
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Figure 4. Left panel: Chandra surface brightness profiles in the selected angular sectors discussed in the text (section ^3.2} and shown in Fig. |3] Black 
circles: sector PA 340° -130°. Red circles: sector P.A. 130°-160°. Green circles: sector PA 160° -250°. Blue circles: sector^0°-340°. The inner edge and 
enhancement of the surface brightness profile are evident in the sector P.A. 130° -160° . A steepening of the sector P.A. 250° -340° profile can also be seen at the 
edge of Chandra field of view, consistent with the XMM feature of Fig[2] Right panel: Zoom over the interesting radial region for the inner edge. The slope of 
the source profile is rapidly changing in the sector P.A. 130° -160° than m the control sector P.A. 340° -130°. 







Table 1 








Parameters from the spectral fits across the edges 








XMM 






CHANDRA 






T 


Zpe 




T 


Zpe 


Radii PA 


X^/dof 


(keV) 


solar 


X^/dof 


(keV) 


solar 


Outer edge 
100- 150^' 250° -340° 


98/105 


^47+0.05 

1 40+0:09 


0.53+0-07 
co+o:o7 

U.D^ 07 


64/70 


1 57+0.05 
i-^/_0 06 


52+0-10 


180-230^' 250° -340° 


55/68 








Control regions 
100-150^' 340°-70° 
180-230^' 340° -70° 


84/122 
73/91 


^59+0.04 

1 42+8:88 
1.48+8-84 
1 58+8:85 

^•-^^-0.06 

I 22+0-04 

1.34+0-01 
1 ,10+0:04 


53+0-07 

0:36+11 

0.49!:| 
0.44!8-87 

0-38!0-07 

41+0-09 
^•^^-0.09 

Q 59+0.11 

0.53+0-07 
050+8:85 

4Q+8:86 
^•"^^-0.06 


58/66 


1 CQ+0.05 
^•-^^-0.06 


57+0-11 
^•-^'-0.10 


100-150^' 160° -250° 
180-230'^ 160° -250° 

Inner edge 
60-90^' 130°-160° 


139/150 
94/113 

48/54 


72/72 
24/21 


1 CQ+0.05 
^•J^_0.06 

1 57+0-08 

2:00^;! 
!:?3i 

1 cQ+8:85 

^•^^-0.06 


55+0-10 
'^•'^''-0.09 

0.45+0.14 

A9+8:i6 

0.53+0-17 
57+0:11 


100-150'' 130°-160°t 

Control regions 

60-90'' 340° -70° 


29/50 
63/81 


39/47 
29/41 


100-150" 340°-70° 


84/122 


58/66 


60-90" 160°-250° 
100-150" 160° -250° 


126/123 
139/150 


60/55 

72/72 



Note. — Results of XMM and Chandra spectral fits, in the radial bins inside and outside the surface brightness edges, as discussed in section Q The first 
and the second columns refers to the radial and angular range of each bins. XMM and Chandra are reported as labeled on the top of the table. No entry in the 
Chandra rows means that the referred bins falls out the available field of view. 
^ The Chandra data have been extracted in a 100° -160° PA region to increase the statistics 



temperature profile along the outer edge is available only in 
the left panel of Fig|6] 

For the inner edge both XMM and Chandra temperature 
profiles show an enhancement of the gas temperature across 
the discontinuity; the same behavior is not observed in the 
control sectors. It is also clear that as expected there is asym- 
metry in the temperature profile, in particular the part of the 
"mushroom head" of emission in the SE leading to the in- 
ner edge is cooler than the surroundings. The temperature 
trend observed by the two instruments is the same, although 
in the outer radial bin the Chandra temperature is 38% higher 
than the XMM one. We interpreted this discrepancy as due 
to a temperatures mixing in the XMM spectrum, induced by 



its larger PSF. To check this hypothesis we compared XMM 
and Chandra temperatures in two wider (P.A: 100° -1600°) 
radial bins at the same radii where we observed the tempera- 
tures discrepancy (100''-150'0 and further away (150''-210'0. 
Coherently with the given interpretation, away from the dis- 
turbed region where the edge is seen, the discrepancy between 
Chandra and XMM temperatureles falls to the 12% level. 

There is no statistically significant change in the temper- 
ature across the outer edge. The quality of the data allows 
us only to measure projected temperatures meaning that the 
true de-projected difference of temperatures across the front 
can be lower. Furthermore when contrasted with the declin- 
ing trend of the control regions at different azimuthal angles. 
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Table 2 

Parameters from the fits of the surface brightness profiles across the edges. 



r (kpc) 



lout 



Pin 
Pout 



M 



Outer edge P.A. 130°-160° 76 0.88±0.09 0.85±0.03 1.28±0.03 1.27±0.10 
Inner edge P. A. 250° -340° 45 0.90±0.05 1.00±0.06 1.65±0.06 1.30±0.13 



0.55 ±0.08 
57+0-10 



Note. — Results from the XMM and Chandra profiles fits, along the outer and the inner edge respectively. The first two columns refer to the internal and 
external power-law index, as discussed in the text, section ^5] We then report the derived density and pressure jumps and the Mach number M. 




Figure 5. XMM residual map, obtained by subtracting to the image of Fig[T] 
the best fit beta model for the data. The center of IC 1 860 is shown by the 
black cross and the position and extent of the surface brightness edges by the 
red arcs. The "mushroom head" discussed in the text is the central region 
with white colour. The image has been processed to remove point sources. 

the picture is consistent with cooler gas inside the edge. A 
simil ar behavior has been see n for the outer e dge in the NGC 
5044 dGastaldello et al.|2QQ9l ) and NGC 5846 ( |Machaceketal.| 
201 1| ) groups. 



5. FITS OF THE SURFACE BRIGHTNESS EDGES 

In order to characterize in a quantitative way the density 
jump across the cold front we modeled the the surface bright- 
ne ss profile acro ss the discontin uity following the a pproach 
of |Rossetti eTaL] p007) and Gas taldello et al.| ( |2QTQ| ). We as- 
sume that the gas density profile is described by a power law 
on either side of the edge, 



n = flout (j;-) 



r<rcf 
r>rcf 



(1) 



where r^/ is the cold front radius and we derive the parameters 
of this model from the fit of the projected surface brightness 
profile, which can be expressed as the integral of the emissiv- 
ity along the line of sight: 



Y. = K 



n\z)A(T,Z)dz, 



(2) 



where K is a constant and we took into account the depen- 
dence of the emissivity of the temperature T and metal abun- 
dance Z. 

We fit the outer part of the surface brightness profile to set 
the external component parameters and we successively de- 



rive the best fit parameters of the innermost part. We applied 
this method to the XMM profile across the outer edge and to 
the Chandra profile across the inner edge. Our results are 
summarized in Table [^ and in figure [^ we plot the surface 
brightness profile and the best fit model respectively for the 
outer (left panel) and inner edge (right panel). 

For the outer edge the agreement between the chosen model 
and the data is good. We find a density jump of riin/riout = 
1.28 ±0.03 that, when combined with the temperature and 
abundance information in the regions inside and outside the 
front reported in Table [T] gives a pressure ratio of Pin /Pout = 
1 .27 ± 0. 10. If this pressure ratio is interpreted as evidence for 
bulk m otion of the cold front then following Vikhlinin et al.] 
( 2001| ) we infer that the front is moving subsonically with a 
Mach number M = 0.55 ± 0.08. 

In the profile of the inner edge a surface brightness bump 
is seen between ~20''and ~50'^ which is not well fitted by 
the adopted two power law model. This discrepancy was ex- 
pected, given the peculiar shape of the source in the sector 
containing the edge: the observed bump corresponds roughly 
to the region where the source is elongated in a narrow lu- 
minous tongue. Except for this peculiar structure, the profile 
inside and outside the edge is well described by the model and 
then internal and external power law indexes are well deter- 
mined. We find a density jump of ntn/riout = 1.65 ± 0.06 and 
using the Chandra temperature and abundance determination 
the pressure ratio is Pin /Pout = 1.30 ±0.13 corresponding to a 
Mach number M = 0.55;^q|4. We will explore the robustness 
of this result with respect to various analysis choices in |6.5[ 

6. SYSTEMATIC ERRORS 

In this section we provide an investigation of the possible 
systematic errors affecting the temperatures and abundances 
across the edges quoted in Table [T] Readers who are uninter- 
ested in these technical details can safely skip to the following 
section. 

6.1. Plasma Codes 

We compared the results obtained usin g the APEC code 
to tho se obtained using the MEKAL code ( [Kaastra & Mewe 
1993[ ) The different implementations of the atomic physics 
and different emission line lists in the plasma codes produce 
no qualitative differences between the fits and no significative 
variations of x^- Variations of the fitted XMM temperatures 
are all within ~5% and well comprised into the statistical er- 
rors; all fitted XMM abundance agrees within ~9%. 

For the Chandra data the use of the different plasma code 
produces a 5% variation of the inner spectral temperature and 
4% variation of the outer one. 

6.2. Bandwidth 
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Figure 6. Left panel :XMM temperature profile along the sectors containing the surface brightness edge. Blue points: outer edge sector, P. A. 250° -340°. Red 
points: inner edge sector, P.A. 130° -160°. Horizontal error bars marks the radial extensions of the bins we used for the spectral analysis.The vertical dashed 
lines marks the radial position of inner and outer surface brightness edges. Right panel: Chandra temperature profile along the sectors containing the surface 
brightness edge. Green points: outer edge sector, P.A. 250° -340°. Magenta points: inner edge sector, P.A. 130° -160°. Horizontal error bars marks the radial 
extensions of the bin for the spectral analysis. The vertical dashed lines marks the radial position of the surface brightness edge as in the left panel: note that the 
more external radial bin of the XMM temperature profile falls outside Chandra field of view and for clarity of the plot the temperature axes do not have the same 
extent in the two profiles. 



^ 0.010 





Figure 7. Left panel :¥\i of the XMM surface brightness profile along the outer edge. The superposed solid line represents the best fit two power law model, as 
discussed in section ^ Right panel: Fit of the Chandra surface brightness profile along the outer edge. The superposed solid line represents the best fit two 
power law model. Note that, as discussed in section ^ the model does not fit the profile in the radial range ~20^^-~50^^ where the source is elongated in a 
narrow luminous tail. 



We explored the sensitivity of our results to the default 
lower limit of the bandpass, £^niin = 0.5 keV. For compar- 
ison we performed spectral fits with £^min = 0.4 keV and 
E^i^ = 0.7 keV. The fitted temperatures are substantially 
unchanged and consistent within the statistical errors, both 
for XMM and for Chandra; an analogous result holds for the 
iron abundances. 



6.3. Variable Nu 

To che ck for possible deviati ons of A^h from the galactic 
value of |Kalberla et al.| ( |2005| ), we perform a spectral fits 
leaving A^h as free parameters. For the XMM spectra the 
fitted A^H values are higher by a factor ~ 3 with respect to the 
Galactic value; the correspondent variations of fitted temper- 
atures are within ^9%; absolute variations of temperatures 



and abundances are however well accounted by the quoted 
statistical errors. 

The parameters of the Chandra spectrum of the region inside 
the inner edge are basically unchanged and the fitted A^h 
agrees with the Galactic value within la. The spectral fits 
of the region outside the inner edge gives instead a higher 
A^H by a factor ~ 4; the correspondent fitted temperature 
(T=1.97;^Q47) however agrees with the quoted one within the 
statistical errors and its lower limit is still consistent with an 
enhancement of the temperature outside the inner edge. 



6.4. Background 

To account for systematic errors in background normal- 
ization, we performed spectral fits allowing it to vary of ± 
5%. Temperatures and abundances are insensitive to this 
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variations in all the considered XMM spectra; variations 
of Chandra temperatures are all within 10%, while iron 
abundances are basically unchanged. 



6.5. Cold front modeling 

The inner cold front is not perfectly modeled by the expres- 
sion given in Eq[T] and there are discrepancies between the 
spectral parameters at either sides of the front between the 
Chandra and XMM values. We explore the sensitivity of the 
derived Mach number with respect to the choice of the fitting 
range of the surface brightness profiles and the adopted tem- 
perature and abundance values. The fiducial measurement is 
derived by using the Chandra spectral values and fitting the 
internal component in the range 1"-AQ" . If we fit the internal 
component in the range !"-?>{)" we obtain a worse fit and a 
Mach number M = 0.35+^ ^ and if we use instead a 60''-90'' 
range to better model the excess, with the problem of deriv- 
ing a very high density in the inner regions, we obtain a Mach 
number M = 0.48;^Q|g. Using the XMM spectral parameters 
instead of the Chandra values we obtain M = 0.67;^q Qg. We 
therefore conclude that our determination in ^ is rather ro- 
bust to possible systematic errors. 

7. OPTICAL ANALYSIS 

In this section we present our optical analysis that includes 
the determination of group membership and substructure de- 
tection. The Digital Sky survey (DSS) image of the central 
region of the group is shown in Figure [8] The purpose of 
the optical analysis is to search for substructures and to use 
them together with the X-ray observations to investigate the 
dynamical state of the IC 1860 group. 

7.1. Sample selection and group membership 

We used the NASA/IPAC extragalactic database (NED) to 
collect galaxies with known recession velocities within 60' 
(1.6 Mpc at the redshift of the source) of the optical position 
of the galaxy IC 1860 which coincides with the peak of the 
diffuse X-ray emission of the group. This corresponds to 1.7 
times the virial radius (rioo) estimate d from hydr ostatic mass 
analysis of the XMM data in Gastaldello et al | (|2()07): this 
is justified by the search for a possible pertur ber following 
the scenario of |Ascasibar & Mark evitch' (1 20061). Th e primary 
sources of these velocities are: 2d FGRS ( Colless et al. 2003), 
the redshift catalogue of 'Dressi er & She ctman' fl988), the 
Southern Sky Redshift Survey (da Costa et al. 1998 ), the red- 
shift survey using the ESO OPTOPUS multifiber instrument 
performed by Stein, ( ,1996 ), the 6dFGS ( .Jones e t al. 2009 ) , and 
the Durham/UKST (DUKST) redshift survey (Ratcliffe et al. 
p^98) . Group members have been selected through the elim- 
mation of background and foreground galaxies along the line 
of sight to IC 1860. An initial rejecti on was performed using 
the "velocity gap" method outlined by |De Propris et al.| ( |2002 ) 
where the galaxies are sorted in redshift space and their ve- 
locity gap, defined for the nth galaxy as Av„ = czn+i -czn, 
calculated. Clusters and groups appears as well populated 
peaks in redshift space which are well separated by velocity 
gaps from the nearest foreground and background galaxies: 
we uses 1000 kms~^ as velocity gap and Fig|9] shows that 
this procedure clearly detects IC 1860 as a peak at z = 0.02 
populated by 81 galaxies. This initial 81 members list for the 
IC1860 group is shown in Table|3j based mainly on the homo- 
geneous 2dF sample which consists of 74 objects. Additional 



objects come from the other catalogues and we checked in 
Appendix [B] the consistency of the various measurements. 

As a further refinement of the membership allocation pro- 
cess we used the slightly mod ified version of the "shift- 



ing gapper", fir st employed by JFadda et al.| ([1996), as de- 
scribed by Owers et al.l (|2009bl). The method^utilizes both 
radial and peculiar velocity information to separate interlop- 
ers from members as a function of group-centric radius. The 
data are binned radially such that each bin contains at least 
20 objects. Within each bin galaxies are sorted by pecu- 
liar velocity with velocity gaps determined in peculiar ve- 
locity. Peculiar velocities were determined by estimating 
the mean group velocity using the biweight location estima- 
tor ( [Beers et"ar. 1990) which was assumed to represent the 
cosmological redshift of the group, Zcos- The peculiar red- 
shift of the galaxy is Zpec = (^gai -^cos)/(l +^cos) and the pe- 
culiar velocity is derived using the special relativistic formula 
Vpec = c((l +Zpec)^ " 1)/((1 +^ ipec)^ + 1), whcre c I S the speed of 
light. The "f pseudosigma" ( [Beers et aL]|1990| ) derived from 
the first and third quartiles of the peculiar velocity distribu- 
tion was used as the fixed gap to separate the group from the 
interlopers, because of its robustness to the presence of inter- 
lopers. The above procedure was iterated until the number of 
members was stable. The results are shown in Fig 10 and the 
rejected galaxies are highlighted in Table [3] 

The final group sample contains 74 members. The value for 
the biweight location estimator of the mean group velocity 
is 6753 ±44 kms~^ which corresponds to Zcos = 0.02252 ± 
0.00015. We used the biweight scale estimator to estimate 
a velocity dispersion of 401 ±46 kms~^. The errors for the 
redshift and velocity dispersion are la and are estimated using 
the jackknife resampling technique. 

In FigflT] we show the velocity dispersion profile derived 
for the data. Bins contain 10 galaxies each and velocity 
dispersions have b een calculated using the gapper algorithm 
( [Beers et aL][1990| ). The profile shows the typical falling be- 
havior with radius expected for a relaxed object. 



7.2. Peculiar Velocity of the Central Galaxy 

We calculated the peculiar velocity of the central galaxy IC 
1860 which is 1 13 ± 50 km s~\ just adding in quadrature the 
errors on the galaxy and group velocities. An offset velocity 
must lie outside a ran ge defined by the cluster v elocity disper- 
sion to be significant ( Gebhardt & Beers|1991| ). We therefore 
calculated th e Z-score for the central galaxy (eq. (7) of Geb- 
[hardt & Be ers 1991 ) which is 0.279 and its 90% confidence 
limits are (0.036, 0.520). The confidence interval has been 
calculated using 10000 bootstrap resamplings and allowing 
for the measurement error of the central galaxy velocity by 
sampling from a Gaussian distribution with standard devia- 
tion set to the reported error as part of the bootstrap. The 
velocity offset is marginally significant at the 2. 3 a level. 



7.3. Substructure Detection 

Given the presence of sloshing cold fronts in the X-ray 
emitting medium of the IC 1860 group it is relevant to look 
for the presence of substructure and to possibly identify the 
merging sub-halo responsible for the onset of sloshing. 

In general the higher the dimensionality of the substructure 
test, the more sensitive it is to substructure; however the sen- 
sitivity of individual diagnostics depends on the line of sight 
relative to the merger axis and therefore no single substructure 
test is the most sensitive in all the situations; it is essential to 
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Figure 8. DSS image of the central region of the IC 1860 group. IC 1860, IC 1859 and ESO 416 - G033 are indicated (see discussion in the text [T33) . IC 
1858 is also indicated (see discussion in Appendix [C). Red labels refer to galaxies rejected by the shifting gapper method (see |7.1) , with numbers refermg to the 
number entry in the catalogue of Table [3] 
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Figure 9. Histogram of all the available redshifts in NED within a 60^ field 
around IC 1860 (galaxies with z > 0.25 are not shown for display purposes). 
The IC1860 group stands out as the redshift peak at z = 0.02. 



apply a full range of statistical tests (e.g.,'Pinkney e t al.|1996| ). 
In this section we apply a battery of statistical methods to the 



Figure 10. Peculiar velocities as a function of the distance from the center 
of the IC 1860 group, illustrating the procedure of refinement of the member- 
ship using the shifting gapper technique (see text). The black circles represent 
galaxies allocated as group members while the red triangles are rejected fore- 
ground and background galaxies lying close to the group in redshift space. 
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Table 3 

IC1860 Group Galaxies 



ID 


Galaxy Identifier 


RA 


DEC 


cz 


cz Error 


Source 


1 


IC 1860 


42.39047 


-31.1891 


6866 


20 


Stein (199^ 
Stein ( lyyg 


2 


2MASX J02493302-3111569 


42.38674 


-31.19923 


7532 


33 


3 


2MASX J02493200-3 110229 


42.38304 


-31.17301 


6865 


64 


4 


2dFGRS S394Z186 
2MASX J02493933-3112169 


42.38183 
42.41373 


-31.16742 
-31.20474 


7525 
7905 


89 
20 


2dF 


5 


|Stein|( 1996| 
zdt 


(5* 


2MASX J02494174-3112139 


42.42371 


-31.20347 


8664 


64 


7 


2dFGRS S394Z183 


42.4035 


-31.2205 


6416 


123 


2dF 


8* 


2dFGRS S394Z176 


42.4285 


-31.19736 


8784 


64 


2dF 


9 


2MASX J02492990-3 109249 


42.37436 


-31.15675 


7945 


64 


2dF 


10 


MCG -05-07-036 


42.34263 


-31.18172 


6835 


64 


2dF 


11 


2dFGRS S467Z174 


42.43247 


-31.15795 


6985 


64 


2dF 


12 


2MASX J02493349-3107159 


42.38908 


-31.12111 


6985 


64 


2dF 


13 


2dFGRS S467Z173 


42.45072 


-31.13592 


6835 


64 


2dF 


14 


2MASXiJ0249572-311036 


42.4886 


-31.17639 


6086 


64 


2dF 


15 


2dFGRS S394Z181 


42.398 


-31.09167 


6266 


64 


2dF 


16 


2MASX J02495627-3 107599 


42.48442 


-31.13339 


7195 


64 


2dF 


17 


IC 1859 


42.26635 


-31.17251 


5966 


64 


2dF 


18 


2dFGRS S467Z168 


42.52021 


-31.13019 


7664 


64 


2dF 


19 


2MFGC 02270 


42.48279 


-31.2875 


6356 


64 


2dF 


20 


IC 1858 


42.28504 


-31.28958 


6056 


89 


2dF 


21 


2dFGRS S394Z189 


42.30958 


-31.06367 


6446 


89 


2dF 


22 


LSBGF416-020 


42.50358 


-31.30597 


6835 


89 


2dF 


23 


2MASX J02485680-3 106431 


42.23634 


-31.11194 


6835 


64 


2dF 


24 


2MASX J0249 1494-3 103060 


42.31196 


-31.05162 


6760 


23 


Stein (1996) 


25 


2MASX J02493567-3 101249 


42.39858 


-31.02364 


6955 


64 


26 


2dFGRS S394Z174 


42.45371 


-31.03258 


7345 


123 


2dF 


27 


2dFGRS S467Z163 


42.56 


-31.29644 


7315 


89 


2dF 


28 


ESQ 416- G 033 


42.40642 


-30.99583 


6595 


123 


2dF 


29 


2dFGRS S466Z088 


42.23542 


-31.38433 


6685 


89 


2dF 


30 


2dFGRS S394Z172 


42.47059 


-30.95476 


7045 


64 


2dF 


31 


2dFGRS S394Z162 


42.68179 


-31.04658 


7045 


123 


2dF 


32 


ESQ 416- G 035 


42.63233 


-31.39464 


6535 


64 


2dF 


33 


2dFGRS S394Z161 


42.68804 


-31.04619 


6625 


89 


2dF 


34 


DUKST 416-040 


42.1205 


-31.37789 


6026 


123 


2dF 


35 


2dFGRSS394Z115 


42.22129 


-30.88644 


6505 


123 


2dF 


36 


2MASX J025 11626-3 119392 


42.81783 


-31.32758 


6559 


45 


6dF 


37* 


ESO416-G025 


42.16992 


-31.53622 


4977 


64 


2dF 


38 


ESQ 416- G 036 


42.64808 


-31.52008 


7015 


64 


2dF 


39 


2MASX J02504193-3130145 


42.67467 


-31.50403 


7135 


64 


2dF 


40 


ESQ 416- G 034 


42.60833 


-31.55828 


7045 


64 


2dF 
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Figure 11. Velocity dispersion profile of the IC 1860 group, where IC 1860 
is assumed to be the center of the group. Bins contain 10 galaxies each and 
Icr jackknife errors are shown. 



7.3.1. One-dimensional Tests for Substructure 

We analyze the velocity distribution to look for departures 
from Gaussianity which can be attributed t o dynamical activ - 
ity. As a first test of Gaussianity, following Hou et alT] ( |2009| ), 
we applied the Anderson-Darling (AD) test as the more re- 
liable and powerful test at detecting real departures from an 
underlying Gaussian distribution. We used the AD test as im- 
plemented in the task ad.test of the package norte st in version 



2.10 of the R statistica l software environmen f^ R Develop 
|ment Core Team|201 \) . When applied to the final group sam- 
ple of 74 objects (the galaxies not highlighted in Table [3]) the 
AD test retu rns a A^ statistic (following the notation of jHou 
et al.||2QQ9[ just A in the ad.test documentation) of 



which corresponds to a /7-value of 0.1059 (computed using 
the modified statistic A^*), therefore consistent with having a 
Gaussian distribution. 

We estimated three shape indicators, the skewness, the kur- 
tosis, and th e scaled tail index, a robust indicator ( |Bird &| 
Beers|1993] ). The values for the sample are 0.429 and 0.386 



for skewness and kurtosis respectively and they show node 
parture from a Gaussian distribution (see Table 2 of [Bird & 



detection of substructure in IC 1 860. 



http ://w w w.r-proj ect . org 
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Table 3 



Continued 



ID 


Galaxy Identifier 


RA 


DEC 


cz 


cz Error 


Source 


41 


ESQ 416- G 027 


42.23883 


-30.79289 


6835 


64 


2dF 




LU5i3 :SU44iF 






oyzD 






4J 


2d^UK5 54t)/zi:)U 


42./086/ 


-31.b30b8 


6026 


89 


2dF 


44 


2d^UK5 54t)/ZiyU 


42.28^96 


-31.68119 


6 /4b 


89 


2dF 


4^^ 


2dtuK:i:^4ooz.wu 


42.U/21 


-Jl.b4b/8 


bblb 


89 


2dF 


46 


2a^uK5 53y4zi:):) 


42.919 


-30.94231 


6446 


123 


2dF 


4/ 


2a^uK:s^4t)/zizi 


42.99bb 


-31.09386 


6386 


64 


2dF 


48 


2d^UK5 54t)/Z134 


42.934b 


-31.48431 


686b 


64 


2dF 


49 


LU^iJ LU144F 


42.83M6 


-30./8181 


/OlO 


31 


6dF 


bO 


2MASA JUZ483i:)U-3 142206 


42.1312b 


-MJUb/b 


683b 


64 


2dF 


bl 


~^MA5A JU24y3334-314:)3iy 


42.3^/2:) 


-31./:)880 


Wbb " 


89 


2dF 


b'l 


ti^U 410- U UZi 


41./8U:J3 


-^LA^M/ 


743:) ~ 


89 


23F 


b6 


2d^UK5 54t)/Z114 


43.1131/ 


-31.32844 


6b0b 


123 


2dF 


b4 


DUK51 41t)-U3U 


43.0/238 


-30.94333 


66bb 


64 


2dF 


bb 


~"2MaSX J 02483040-3 1473B5~ 


42.126/b 


-31./93b6 


63b6 


64 


23F 


bb 


2d^UK5 54t)t)ZU8 / 


42.244/9 


-31.83344 


69bb 


89 


2dF 


b'/ 


Zai^UK^S ^4t)t)Z14U 


41.79083 


-31.606/8 


6b9b 


89 


2dF 


b8 


2MA5A1 JUZ:)U100-3U3222 


42.b69b 


-30.b39b 


/0/b 


64 


2dF 


b9 


2d^UK5 54t)/Z13:) 


42.91488 


-3l./38b3 


bWT~ 


123 


2dF 


bU-^ 


ZdtuK^ :sj9Jzuj/ 


4LbH4/^ 


-JU. /^SJO 


4/U/~ 


123 


2dt 


6i 


DUK5 1 41b-U43 


41.638/b 


-31.M614 


/34b 


89 


2dF 


62 


h5U 410- U U4U 


43.11304 


-30.//6b3 


6//b 


89 


2dF 


63 


2d^UK5 540 /Z13/ 


42.8M 


-31.82 


6//b 


89 


2dF 


64 


2d^UK:S 53y4ZU41 


42.9261/ 


-30.b9/b 


6 /4b 


89 


2dF 


bb 


2d^UK5 5400Z1U1 


42.14604 


-31.911b3 


686b 


89 


2dF 


66 


Li^^a 5U44UF 


41.86242 


-31.81b89 


6838 


IbO 


DUKSr 


6/ 


AMUZ:)U-31U 


43.22633 


-30.8b /42 


6 /4b 


64 


2dF 


6^=^ 


IdtUK^ :sJ9Jzuz:) 


4L8W/^ 


-MbUUJb 


b276 


89 


2dt 


69 


2d^uK5 53y:)zuz / 


43.0041/ 


-30.b664/ 


69bb 


89 


23F 


/U 


2(l^UK5 540/Z130 


42.8 /2b8 


-31.92b83 


6146 


123 


2dF 


7i 


2d^UK5 53y:)ZU3 / 


43.33496 


-30.914b6 


6/ lb 


89 


2dF 


/2 


b5U 410- u u^y 


43.09304 


-31.81/28 


662b 


64 


23F 


/3 


2d^UKb 53y4ZUn 


43.2806/ 


-30.//b33 


6416 


64 


2dF 


74 


2d^UK:S ^i^iZU^U 


41.80092 


""^U.4/830 


08U3 


89 


2dF 


lb 


2d^uK5 53y:)ZU3:) 


43.41 //9 


~^ 1.10042 


0380" 


64 


2dF 


lb 


hbU410-UU41 


43.3/23/ 


~^3U.80U30 


0410~ 


64 


2dF 


//r 


latuK^i :i400Aizz 


41.^8^21 


-S2mii9 


494/ 


b4 


2dt 


/8 


L5tJU^410-U18 


43.264 


-31./b99/ 


6296 


123 


2dF 


/9 


2d^UK5 5400Z183 


41.43^1/ 


-3i.oy:)z:) 


0!)^:) ~ 


89 


2dF 


8U 


2d^UK5 5400Z10/ 


41.^^2Zb 


-31.8b809 


0446~ 


123 


2dF 


8i 


2d^UK5 5400Z100 


41.bb6b8 


-31.88614 


698b 


89 


2dF 


1 



Note. — C olumns are as follows: (1) Galaxy Identifier (2) Right ascension (J2000), (3) Declination (J2000), (4) Recession Velocity, (5) Error on the Recessio 
Velocity (^Source of the measurement 
Ihe galaxies with entry in italics and marked with an asterisk are the galaxies rejected trom memebership by the shitting gapper method. 



Beers 1 199? , for 75 data-points). The scaled tail index is 1.14 
again consistent with a Gaussian distribution. 

We investigated the presence of significant gaps in the ve- 
locity distrib ution following the weighted gap analysis of 
Beers et aL] p991) looking for normalized gaps larger than 
2.25 since in random draws of a Gaussian distribution they 
arise at most in about 3% of the cases, independent of the 
sample size. We did not find any significant gap in the sam- 
ple. 

To detect subsets in the velocity distribution we als o ap 



plied to the data the Kaye's mixture model (KMM) test ( [Ash- 
man et al^|1994 ). The KMM algorithm fits a user-specified 



number ol Gaussian distributions to a dataset and assesses the 
improvement of that fit over a single Gaussian. In addition, it 
provides the an assignment of objects into groups. The KMM 
test is more appropriate in situations where theoretical and 
empirical arguments indicate that a Gaussian model is rea- 
sonable, as for the line of sight velocities in a dynamically 
relaxed cluster. We did not find again any significant partition 
of the sample. 
To summarize the results of the ID tests, we found no ev- 



idence of substructure in the velocity distribution in the fi- 
nal group sample. Note that the same battery of tests applied 
to the unrefined sample would detect a significant departure 
from the Gaussianity. Thus, the shifting gapper refinement 
of the membership allocation has been effective in remov- 
ing small foreground and background groups; in particular the 
ce ntral group of high ve locity interlopers was already noticed 
by [Burgett eFaTl i l^O?! ). 



7.3.2. Three-dimensional Tests for Substructure 

The existence of correlations between positions and veloci- 
ties of cluster galaxies is a signature of true substructures. As 
a first approach we adopt the Dressler-Shectman A statistic 
( [Dressier & Shectma n 1988 ) which tests for differences in the 
local mean and dispersion compared to the global mean a nd 
dispersion and it is recommended by'Pinkn ey et al.[ ( [l996| ) as 
the most sensitive 3D test. Calculation of the A statistic in- 
volves the summation of local velocity anisotropy, 5, for each 
galaxy in the group, defined as 
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Figure 12. Results for the k substructure test. Each dot represents a galaxy 
in the IC 1 860 group and the size of the circle around the dot gives an indica- 
tion of the difference of the local velocity distribution compared to the global 
group velocity distribution (the size is scaled by hit where hit is given by Equa- 
tion Bj. Clustering of large circles would indicate significant departures and 
the liKely presence of a substructure. 



S^ = 
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ocal ^groupj v<^local ^) J? 



(3) 



where A^nn is the number of nearest neighbours over which 
the local recession velocity (viocai) and velocity dispersion 
(c^iocai) ar e calculated . Here we adopt A^nn = V^, following 
[Pinkney e t al.| ( |199 6|). The significance of A is estimated 
using 10000 Montecarlo simulations randomly shuffling the 
galaxy velocities. We find suggestions of substructure with 
the DS test at the marginal level 96.7% level. A similar re- 
sult (95.8%) is obtained when using the k test of Colless &| 
Ipunn (1996). The A>:-test searches for local departures from 
the global velocity distribution around each cluster member 
galaxy by using the Kolmogorov-Smimov (KS) test which 
determines the likelihood that the velocity distribution of the 
V^ nearest neighbors around the galaxy of interest and the 
global cluster velocity distribution are drawn from the same 
parent distribution. Calculation of the the k, statistic involves 
the summation of the negative log likelihood for each galaxy 
defined as 



f^i = -logPKsiD > Dobs.i), 



(4) 



over the entire sample. By inspecting the bubble plot of 
Fig[T2]the galaxies contributing to the signal of the n statistic 
are me galaxies with high velocities in the center of the group 
(including the galaxy IC 1860 itself) and a subclump in the 
north-east outskirts of the group (the galaxies with ID # 75, 
ID # 76, ID # 73, ID # 81, ID # 54, and ID # 56). Remov- 
ing these galaxies and re-running the test results in a A value 
being significant at only the 71% level. If on the contrary we 
remove the galaxies identified in the center with high values 
of 5 and we re-run the test we obtain a A value being signif- 
icant at the 81% level. It is important to bear in mind that 
the A stati stic is particularly insensitive to superimposed sub- 
structures ( [Dressier & Shectman||l988t |Pinkney et aL||1996| ) 
and that it has the highest false positive detecti on of substruc- 
tures due to velocity gradients across the group ( [Pinkney et al. 



[T9961 ). 

We tested further the presence of these two candidates sub- 
groups by applying the full 3D-KMM method. The 3D-KMM 
test, starting from the above mentioned seeds, can allocate 
galaxies in subclusters. Using the NE subgroup and the rest 
of the galaxies as initial two-groups partition, we did not find 
any significant allocation. 

On the contrary a two-group partition with as initial seeds 
the inner subgroup and the rest of the members we find that a 
partition of 60 and 14 galaxies is a more accurate description 
of the 3D distribution than a single Gaussian at more than the 
99.999% level. 

To summarize the results of the 3D tests we found marginal 
evidence of substructure related to a possible subgroup in the 
NE outskirts of the group and a central structure with rela- 
tively high peculiar velocities. 

7.3.3. Candidate Perturbers 

The possible substructures identified by the statistical tests 
applied above do not seem likely to be the perturbers. The 
central substructure comprises the central galaxy, IC 1860, 
itself and the north-east substructure is a loose aggregate of 
very faint and small galaxies. 

A visual inspection (see FiglSl) revealed two possible can- 
didates: IC 1859 and ESO 416^0033. IC 1859 is a Seyfert 
2 galaxy with a peculiar optical morphology showing a pair 
of spurs in his easte rn spir al arm; it was already classified 
as peculiar in Kelm et al.[ ( p^998) . It is at a projected dis- 
tance of 172 kpc from IC 1860 (corresponding to 0.18 ryir, 
|Gastaldello et aL][2007 ) and it has a velocity difference of 
-900 ±67 kms-' with respect to IC 1860 and of -787 ± 78 
kms~^ with respect to the group mean velocity. ESO 416- 
G033 has a mildly disturbed optical morphology, in particular 
if compared to IC 1859, and it is at a projected distance of 
310 kpc (0.33 Tvir) from IC 1860 and it has a velocity differ- 
ence of -271 ±124 kms~^ with respect to IC 1860 and of 
-158 ± 130 kms~^ with respect to the group mean velocity. 
A third galaxy SW of IC 1858, the SO IC 1858, is a striking 
narrow-angle tail radio galaxy (NAT) as shown by the GMRT 
data and discussed in Appendix O Thanks also to the radio 
tail constraining the motion in tne plane of the sky we are 
able to rule out this galaxy as the perturber (see ^. 

8. RADIO ANALYSIS 



A r adio image of IC1860 has been presented by [Dunn 
[et al.[ (2010), who analyzed Very Large Array (yiA) A- 
configuration data at 1.4 GHz with angular resolution of ~ 3". 
A point source with two faint, small-scale extensions was 
found associated with the galaxy, but no correlation between 
these radio features and the X-ray emission was observed. 

In search for possible extended radio emission, undetected 
in the current radio images, we analyzed low radio-frequency 
observations of IC1860 extracted from the Giant Metrewave 
Radio Telescope (GMRT) archive (project 06EFA01). The 
source was observed at 325 MHz in September 2004 for ap- 
proximately 8 hours on source, using only the upper side band 
(USB) with 16 MHz bandwidth. We calibrated and reduced 
the data set using the NRAO Astronomical Image Processing 
Sy stem package (AIPS), fol lowing the procedure described 
in [Giacintucci et al.[ ( [2008[ ). We applied phase-only self- 
calibration to reduce residual phase variati ons and improve 
the quality of the final images. The scale of Baar s etal.[ ( [1977 ) 
was adopted for the flux density calibration. Residual ampli- 
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Figure 13. GMRT radio contours at 325 MHz overlaid on the map of XMM 
residuals (FigjSJ. The restoring beam is 30^^ x 29^^ in p. a. 6° (white ellipse). 
The rms noise level is 0.6 mJy beam~^ . Contours are spaced by a factor of 2 
starting from 2 mJy beam~^ . The emission which can be seen in the South is 
coincident with an X-ray point source. The inset shows the 325 MHZ image 
at higher resolution (18^^ x 14^^ in p.a. -24 °; black ellipse), overlaid on the 
Chandra image. The rms noise level is 0.5 mJy beam~^ . Contours start at 1.5 
mJy beam~^ and then scale by a factor of 2. No negative levels, corresponding 
to -3a, are present in the portion of the image shown. 



{■ 



tude errors are within 8% (e.g., Chandra et al.|2Q04| ). 

Fig 13 shows the GMRT contours at 325 MHz overlaid on 
the map of XMM residuals shown in Fig|5] To highlight the 
extended emission, the radio image has been produced with a 
resolution of ^ 30'^ obtained by applying a taper to the uv- 
data through the parameters ROBUST and UVTAPER in the 
task IMAGR. The radio image of the central region is also 
shown at higher resolution (18"), overlaid on the Chandra 
image. The image shows an unresolved central component, 
spatially coincident with the emissio n detected at h igher fre- 
quency and higher resolution by Du nn et al. ( 2Q1Q| ). Its flux 
density at 325 MHz is 60.3 ±4.8 mJy, as measured on the 
full resolution image (14^^ x 9"', not shown here), and a total 
flux density of 13.9 mJy is reported at 1.4 GHz by [Dunn et al. 
2010), implying a spectral inde^xQJa =1.0. 

In addition to the central point source, the GMRT image 
reveals an extended feature with very low surface brightness 
which originates at the IC1860 center and extends ~ 100^' (45 
kpc) from the center toward south-east. This radio emission 
is contained within the inner cold front and traces the "mush- 
roo m he ad" feature interpreted as the tip of the sloshing spiral 
in |3.3[ A flux density of '^22 mJy is associated with this 
component. A second patch of faint radio emission is visible 
~ 2' south-west of the central source, with a flux density of 
~20 mJy. Neither of these e xtended features are visible in the 
image by Dun n et al. ( |2010| ), due to the lack of short spacings 
in the uv coverage of their data set (the largest angular scale 
that can be imaged with the VLA in A configuration at 1.4 
GHz is - 3 8^0. We thus inspec ted the NRAO VLA Sky Sur- 
vey (NVSS, [Condonetal.|1998l ) image at 1.4 GHz, but found 
no evidence of significant extended emission over the two ex- 
tended structures detected by the GMRT. Given the sensitivity 
level of the NVSS image (la = 0.45 mJy beam"^ one beam 
corresponding to 45^0. we placed a lower limit of a > 1.9 
for the spectral index of the SE extension and a> 1.8 for the 
SW emission. Both spectral index limits are very steep and 

^^ We adopt the convention Su oc u~^ for the synchrotron spectrum, where 
Siy is the flux density at the frequency v. 



suggest that the extended emission originates from aged radio 
plasma, perhaps associated with a past episode of activity of 
the central radio galaxy. 

9. COMPARISON WITH SIMULATIONS 

The analysis of the available Chandra and XMM data for 
the IC 1860 group has revealed a pair of surface brightness 
discontinuities at 45 kpc and 67 kpc that the spectral analysis 
indicated to be cold fronts. We also detected a spiral- shaped 
excess in the core with the outer edges of this spiral being 
the cold fronts. These features are all indicators of ongoing 
slos hing in the scenario put forward to explain these phen om- 



ena CAscasibar & M arkevitc h|2006 
Following Roediger et al. 



i 



IRoediger et al.pon) . 
,011|l2012b| ), we use the slosh- 



ing signatures to derive the merger history. We perform a 
qualitative comparison to the sloshing simulations Roediger 
et al. (2011) performed for the Virgo cluster. The sloshing 
structure appears clearest in brightness residual maps and we 
show such maps derived from the simulations for different 
lines of sight (LOSs) in Figure [14] The top row shows the 
residual map for the LOS perpenmcular to the orbital plane. 
In the left-hand- side column the cluster is rotated around the 
horizontal axis, in the right-hand side column around the ver- 
tical axis. When the LOS is perpendicular to the orbital plane, 
the residual map shows the well-known spiral morphology. 
This morphology can be clearly recognized for inclinations of 
the orbital plane up to 45° away from this configuration. Be- 
yond 70° the spiral morphology becomes weak. The residual 
map of the IC 1 860 displays a spiral-like morphology, there- 
fore we conclude that our LOS is within 70° of the face-on 
view. We note that only the rotation around the horizontal 
axis in Fig. [14] (i.e. the cases in the left-hand- side column) 
leads to a substantial radial component for the orbit of the 
perturber, i.e. a LOS velocity component. The configurations 
in the right-hand- side column lead to a progressively smaller 
angle between the LOS and the orbital plane, but the LOS 
is always perpendicular to the orbit itself. A characteristic 
feature of all cases in the left-hand-side column and the top 
panels in the right-hand- side column is the pronounced large- 
scale asymmetry, with a clear brightness excess in the right, 
i.e. at the side of the cluster center opposite the perturber's 
pericentre. In all these cases, also in projection the pericentre 
of the perturber orbit is left of the cluster center. Only in the 
2 or 3 bottom panels of the right-hand side column, where the 
pericentre is seen more and more behind/in front of the clus- 
ter center, this large-scale asymmetry is weak or absent. In 
IC 1860, the brightness outside the radius of the western cold 

front does not seem stro ngly asymmetrical. 

The simulations (e.g.,'Asca sibar & Markevitch|2006| ) show 
that following the spiral structure from the outside inward 
indicates the direction of the infall of the perturbing object. 
Given that the direction of the spiral inflow is clockwise the 
perturber also has a clockwise orbit. We find two configu- 
rations of the orbit projected onto the sky fit the observation 
best. 

E-N-W: The perturber approaches approximately from the 
east, passes the group center in the North and departs 
approximately towards the West (case E-N-W). The 
IC 1860 observation matches this configuration if IC 
1860's western cold front corresponds to the excess re- 
gion to the left in t he top panels in Fig. [T4|( what was 
called cold fan by Roe diger et al.| 20 11), and IC 1860's 
south-eastern cold front to the outermost cold front to 
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Figure 14. Brightness residuals with respect to azimuthally averaged brightness. Original image divided by averaged. Only the inner 500 kpc of the cluster have 
been used to create these images. This slightly affects the quantitative values of the brightness residuals at periphery of the images, but they are qualitatively 
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the right of the core in the to p panels in Fig. 141 wh ich 
was called major cold front in |Roediger et al. POl 1| ). In 
the Virgo simulations, this cold front moved outwards 
with a velocity of '^55 kpc per Gyr. Assuming that 
in the IC 1860 group the cold font moves at compa- 
rable velocity, the group-centric distance of 46 kpc of 
this cold front indicates that the pericenter passage of 
the perturber took place about 1 Gyr ago. This config- 
uration predicts a large-scale brightness excess to the 
south. 

S-E-N: The perturber approaches approximately from the 
south, passes the group center in the east and departs 
approximately towards the north (case S-E-N). The 
IC 1860 observation matches this configuration if IC 
1860's western cold front corresponds to the major cold 
front (outermost to the right in top panels of Fig. [14]), 
and the SE front to the secondary front just leffof 
the simulated cluster core. A close inspection of IC 
1860's brightness residual map suggests that the spiral- 
shaped brightness excess may continue beyond the ex- 
cess just inside the cold front in the west, it may con- 
tinue outwards anti-clockwise over north to the north- 
east, maybe even east. The faint outer eastern end of the 
brightness excess spiral would correspond to the bright- 
ness excess fan in the left of the simulated images. If 
the western, now identified as major, cold front expands 
at the same velocity as the major cold front in the Virgo 
simulations, its group-centric distance of 76 kpc indi- 
cates that the pericenter passage of the perturber took 
place about 1.5 Gyr ago. This configuration predicts a 
large-scale brightness excess to the west. 

Correspondingly, orbits W-S-E and N-W-S do not match the 
structure of the brightness residual map of the ICM 1860 
group well, and we regard them as unlikely. Consequently, 
we expect the perturber in the approximately in the north or 
the west or in between, and the pericentre passage has taken 
place 1 to 1.5 Gyr ago. The discussion above is based on 
the LOS being perpendicular to the orbital plane, but is valid 
to configurations at least 45° away from this. Only beyond 
70° the spiral morphology becomes weak. 

The group gas temperature is approximately 1.5 keV This 
corresponds to sound speed of approx. 600kms~^. If we as- 
sume an average relative velocity of the perturber between 
Mach 1.5 and Mach 2, within 1 Gyr it can move around 1 
Mpc away from the group center in real space. 

The optical image in Fig. [8] shows three giant galaxies as 
potential perturber candidates in the IC 1860 group. There 
are two big distorted spirals in N and W, which fit either orbit 
configurations E-N-W or S-E-N. IC 1859 in the W has strong 
radial velocity (~ 800kms~^) with respect to the group, and 
is quite close to IC1860 in projection(176 kpc). In the north, 
ESO 416 has only mild radial velocity (160 and 270 with 
respect to IC1860 or group mean, respectively), and is also 
slightly further away (310 kpc). If we assume both galaxies 
are about 1 Mpc from the group center in real space, it means 
the orbit of ESO 416 makes an angle of 70° with sky, the orbit 
of IC 1859 an angle of 80° to achieve correct projected dis- 
tance. According to the simulated brightness residual maps, 
this is just within the range of inclinations where the central 
brightness excess takes on a clear spiral structure. If the galax- 
ies were passing the group at high velocity, their real-space 
distance and thus the inclination angle between their orbit and 



the plane of the sky may be somewhat smaller. On the basis 
of the current X-ray data, either of these galaxies may have 
been the perturber. A clear detection of the large-scale bright- 
ness asymmetry would be helpful to distinguish both cases. 
The SO galaxy in the south is a NAT with the tail pointing 
southwards, indicating a projected motion towards the group 
center (see Appendix O. Therefore, both, its position and its 
possible orbit disfavorthis galaxy as the perturber. 

10. DISCUSSION 

10.1. The sloshing scenario at the group scale: cold fronts, 
peculiar velocities ofBCGs and perturbers 

The evidence presented in this paper indicates the presence 
of cold fronts in the galaxy group IC 1860 associated with a 
spiral excess in surface brightness in its core. These features 
fit nicely in the sloshing scenario put forward by |Ascasibar| 
& Markevitch (2006). Furthermore, as qualitatively shown 
in ^ the observed geometry of the sloshing features (cold 
fronts and surface brightness excess) can be reproduced by 
simulations and can constrain the recent merger history , as 
done quantitatively for the Virgo cluster and A496 ( Roediger| 
et al. 20111 |2012b|). Simulations tailored specifically to IC 



1 860 and in general to the groups showing sloshing will allow 
another step forward in interpreting the data. 

Another tell-tale sign of the off-axis merger with a satellite 
is the presence of a peculiar velocity of the central galaxy, as 
also suggested by |Ascasibar & Markevitch ( |2006| ). The ellip- 
tical galaxy IC 1 860, the BCG of the group, has a peculiar ve- 
locity of 1 13 km s~^ with respect to the group velocity distri- 
bution, a difference significant at the 2.3cr level. If we identify 
the BCG position with the bottom of the potential well (i.e. 
the dark matter peak), the order of magnitude of the measured 



velocity is consistent with Fig. 4 of |Ascasibar & Markevitch 
(2006) where one can roughly estimate a DM peak velocity 
of ~100kpc/Gyr i.e ~ lOOkms"^ It would be interesting to 
estimate this parameter from simulations also for a group size 
main halo and for different merger geometries. The presence 
of BCG peculiar velocities is another piece of evidence in fa- 
vor of the sloshing scenario: for example an alternative to dis- 
placing the ICM by gravitational interaction with a subclus- 
ter is t h e passage of a shock as proposed by |Churazov et al. 
( |2003| ). [Roediger et al.| ( |2011| ) showed in the case of Virgo 
that the bow shock of a fast moving galaxy causes sloshing 
and leads to very similar cold front structures as observed. 
However in this scenario no BCG peculiar velocity should 
be observed. The same minor mergers which are a common 
phenomena in the growth of galaxy groups and clusters are 
therefore the origin of both the sloshing features in the X-ray 
emitting gas and the peculiar velocities of the central galax- 
ies, now consistently found in massiv e halos (> 10^^ Mq) as 
for example in van d en Bosch et al.| p005 ). van den Bosch 



et al. (2005 ) in-fact suggest as the preferred explanation for 
the latter finding the scenario of a non-relaxed halo, meaning 
that the BCG is located at rest with respect to the minimum 
of the dark matter potential, but the dark matter mass distribu- 
tion is not fully relaxed and the potential minimum does not 
coincide wi th the bar y centre. This is the same scenario put 
forward by A scasibar & Mar kevitch ([2006) (and envisaged by 
Miralda-Escude (1995)): therefore botE" sloshing and peculiar 
velocities of the BCGs are a natural outcome of a hierarchical 
picture of structure formation. 

The key ingredient in the sloshing scenario is the presence 
of a perturber. The identification of the perturber is straight- 
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forwa rd in case of double gaseous components as in Abell 
1644 ( [Johnson et ar|[20T0l) and in the groups NGC 7618 - 
UGC 12491 ( [Roediger et al.|2()T2a| ); it is also identified in a 
convincing manner by means of optical substructure investi- 
gations (e.g. Owers et al. 2011a). However in some cases, as 
for example A496, notwithstanding many observed X-ray fea- 
tures, in particular the key presence of a large scale asymme- 
try, consistent with a dedicated hydrodynamical simulation, 
it was not possible to identify the perturber. It may be im- 
possible to ever identify it because tidal forces can disperse 
the galaxies in the subcluster which might not be recogniz- 
able as a compact structure anymore (Roediger et al. 2012b). 
At the group scale we might have a cleaner test case in this 
respect for the sloshing scenario than in more massive clus- 
ters: in a group environment the halo of a single moderately 
massive spiral galaxy can be the responsible perturber. In ^ 
we explored two possible scenarios consistent with the iden- 
tification of two big distorted spirals, ESO 416 and IC 1859. 
Given the optical appearance of tidal features pointing in the 
projected direction of the central galaxy IC 1860, a preference 
for the most likely candidate can be given to IC 1859. 

This interpretation is strengthened by the analogy with two 
other groups reported to host sloshing cold fronts, NGC 5044 
and NGC 5846. These two groups also have an optically dis- 
turbed spiral galaxy which can be identified as the possible 
perturber. For the NGC 5044 group we highlighted in |A.1| 
a key feature not reported in previous analysis: the presence 
of a large scale asymmetry in terms of a brightness excess in 
the east (see FigfTS]). From their sloshing simulations for the 
Virgo cluster, Roediger et al. (201 1) predicted the existence of 
a large scale asymmetry. In Virgo the asymmetry could not be 
confirmed observationally, due to incomplete data coverage 
and to the fact that Virgo is still a very dynamic cluster where 
sloshing-related asymmetry may be difficult to be disentan- 
gled form other perturbations. In the cluster A 496 for the 
first time the connection between gas sloshing and the pres- 
ence of a large scale asymmetry was confirmed: the asymme- 
try wa s found in simulations a nd observations at the same po- 
sition ( [Roediger et al.|2 012b ). The galaxy group NGC 5044 
is therefore another example which confirms the picture of 
gas sloshing. Given the presence of almost concentric arcs on 
opposite sides of the cluster core (the identified cold fronts) 
and the orientation of the large scale asymmetry, it will be 
possible to constrain quantitatively the merger scenario. Al- 
ready qualitatively the LOS can be constrained to be almost 
parallel to the orbital plane of the perturber and consistent 
with a position East of the center of the group. All these con- 
straints are compatible with the identification of the perturber 
as the tidally disturbed s piral galaxy NGC 50 54 (see Fig[T8]), 
as already suggested by [David et al. ( 2009[ ). For the NGC 
5846 group the characterization of the cold fronts, the possible 
merger geometry and the identification with the optically (and 
HI, [Higdon et al. 1998) disturbed ringed spiral galaxy NGC 
5850 ( see Fi g[20[) have already been discussed in Machacek] 
eTaL] ( [20TT] ). In pl2|we reported the XMM MOS surface 
brightness residual map which is c onsistent with the Chan- 
dra one of [Machacek et al.[ ( [201 1| ) and reported a possible 
evidence of a peculiar velocity of the BCG NGC 5846 with 
respect to the mean group velocity, as observed in other cases 
of sloshing systems. The analysis and interpretation may be 
challenging even if more data are collected given the possi- 
ble identification of the NGC 5846 group as a subgroup of a 
more extended system comprising also the NGC 5813 group 
(Mahdavi et al. 2005| ). 



X-ray sloshing features are witnessing the interaction of the 
perturber with the group halo. If the perturber as in the case of 
the three groups discussed above is a optically disturbed spiral 
galaxy, further constraints might in principle be obtained by 
comparisons of the optical distorsion of the spiral arms pro- 
duced in simulations following the evolutio n of disc galax- 
ies within the global tidal field of the group ( [Villalobos et al. 
,2012, ). 

10.2. The low frequency radio emission in IC 1860 

Mini radio halos, faint diffuse radio emission with radius 
of 100-300 kpc and a steep spectrum (a ~ 1.3), are present 
in a number of relaxed cool-core clusters. The central active 
galaxies of these clusters are not directly powering the diffuse 
radio emission as the radiative time scale of the electrons re- 
sponsible of the observed emission is much shorter than the 
time required for these electrons to diffuse across the cool- 
ing region (Brunetti 2003). Reacceleration of pre-existing, 
low-energy electrons electrons in the ICM by turbulence in 
the core region is a possible mechanism responsible for the 
radio emission. Sloshing m otions can produce significant tur- 
bulence in the cluster core ( [Fujita et "aL][2004[ [ZuHone et al. 
|2013| ) and this solution is strengthened by the spatial corre- 
lation between radio mini-halo emission and cold fronts in 



some clusters hosting mini-radio halos ( [Mazzotta & Giacin- 
tucci 2 008 ). 

An ongoing statistical study of the properties of clusters 
with and without a minihalo shows that minihalos tend to oc- 



cur on ly in very massive, cool-core clusters ( Gia cintucci et al. 
2013[ ). Indeed, no clear detection of a minihalo at the group 



scale has been reported so far. In IC 1860 we have indication 
of a spatial coincidence between a cold front and extended 
radio emission, similar to what is seen in some minihalo clus- 
ters. The GMRT radio emission at 325 MHz is much more 



extended than the 1 .4 GHz emission presented in Dunn et al. 
(2010) which shows a clear point source with faint exten- 
sions m the north-east and south-west directions with no clear 
correlation with the X-ray extension towards the south-east. 
Instead one component of the 325 MHz emission traces the 
bright tip of the excess spiral feature and is confined within 
the inner cold front, very much like the spatial coincidence 
found in clusters with mini radio-halos. Another component 
is more detached lying in projection at the edge of the spiral. 
This latter feature resemb les the detached c omponen ts found 
at 235 MHz in NGC 5044 ( [Giacintucci etaL|201 1,. David et al.| 
[2009] ). 

However the extension (only 45 kpc) and the ultra- steep spec- 
tral index (a > 1 .9) of the radio emission can be explained 
without the need for in situ re- acceleration as for radio ha- 
los and radio mini-halos in galaxy clusters. Indeed, as dis- 
cussed in Jaffe (1977), a possibility for the transport of the 
electrons responsible for the radio emission is l arge scale mo- 
tions o f the gas carrying the electrons along. [ZuHone et al. 
( 2013| showed how relativistic electrons preferentially lo- 
cated m AGN-blown bubbles are redistributed by the sloshing 
motions so that most of them end up within the spiral shape 
traced out by the cold fronts. Observations at different radio 
frequencies and simulations tailored for the groups mass scale 
are needed to investigate the possibility that re-acceleration 
might take place. 

10.3. Reconstructing the mass distribution in a sloshing 

group 
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Since the introduction of the sloshing paradigm in the 
cluster Abell 1795, it was pointed out the the presence of 
cold front discontinuities in the ICM density and temperature 
profiles can caused flaw ed results in cluster mass estimates 
([Markevit ch et al.||2QQl| ). Unphysical discontinuities arise at 
the cold front in the derived cluster mass profile assuming hy- 
drostatic equilibrium. This is true for most of the published 
cold front profiles and this fact seems consistent with tangen- 
tial flows associated wit h sloshing con tributin g an additional 
centrifugal acceleration ( [Keshet et al.||2010p . [Roediger et al.' 
(20 12b) address the issue to what extent the underlying mass 
distribution can be inferred from fully azimuthally averaged 
temperature and density profiles with the presence of cold 
fronts in their simulations for the case of the Abell 496 cluster. 
They found that sloshing modifies the azimuthally averaged 
profiles by less than 10% and concluded that given the quality 
of the current available data, the distortion of the profiles due 
to sloshing does not seem to introduce an observable differ- 
ence and will not be a major error source for clust er mass pro- 
file. This is true also at the group mass scale: in [Gastaldello 
et al.| ( |2007| ) this issue was addressed directly in IC 1860 by 



comparing the mass profile reconstructed in two sectors with 
PA 15°-95° and 195°-15° finding results on the derived total 
mass (and NFW concentration parameter) within their 2cr er- 
rors. Therefore the low-level disturbances due to sloshing do 
not indicate a significant violation of hydrostatic equilibrium. 

11. SUMMARY AND CONCLUSIONS 

We have presented results of a combined analysis of the cur- 
rently available XMM and Chandra X-ray data, optical spec- 
troscopy and radio GMRT for the bright nearby galaxy group 
IC1860. The results can be summarized as follows: 

1 . the X-ray imaging and spectroscopy reveal the presence 
of a set of two cold fronts at 45 and 76 kpc and of a spi- 
ral of excess surface brightness. Those features are con- 
sistent with the sloshing scenario put forward to explain 
these features, as further supported by the presence of 
a peculiar velocity of the central galaxy IC 1 860 and 
the identification of a possible perturber in the optically 
disturbed spiral galaxy IC 1859; 

2. the analogy of the IC 1860 group with two other X-ray 
groups showing cold fronts, NGC 5044 and NGC 5846 
is tantalizing in terms of the identification of the possi- 
ble perturber with an optically disturbed spiral galaxy 
and peculiar velocity of the bright central galaxy. We 
reported further evidence for the sloshing scenario in 
those two systems. We detected a large scale surface 
brightness excess i n NGC 5044 which strengthen ear- 
lier suggestions by |David et aL] ( |2009| ) of the identifi- 
cation of the perturber with the disturbed spiral NGC 
5054. We detected hint of a peculiar velocity of the 
central galaxy NGC 5846. 

3. The presence of X-ray sloshing features, peculiar ve- 
locity of the BCG and optically disturbed features in 
the perturbing spiral galaxies are signatures of the in- 
teraction of the group tidal field with a disk galaxy and 
indication that also this process is actively participating 
in the transformation of late type galaxies in the group 
environment. 

4. The low frequency GMRT radio data show extended 
emission with a component confined by the inner cold 



front as seen in some clusters with mini-halos. Unlike 
minihalos that require in situ particle reacceleration, the 
extended emission in IC 1860 is consistent with aged 
radio plasma redistributed by the sloshing motions. 

The case studies presented in this paper show the rich phe- 
nomenology and the constraining power that X-ray, optical 
and radio data can provide in assessing the dynamical state 
and past merging history of galaxy groups and their galax- 
ies. We have provided the strongest evidence to date that cold 
fronts in groups exist and have similar properties to those in 
massive clusters. This is critical confirmation of the currently 
favored model of cold- fronts formation in relaxed objects. 
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A. THE NGC 5044 AND NGC 5846 GROUPS 

In this section we report additional X-ray and optical analysis for the NGC 5044 and NGC 5846 groups, which show sloshing 
cold fronts. The analogy of IC 1860 group with these systems will reinforce the suggestion of an optically disturbed spiral galaxy 
as the candidate perturber and provide further evidence of the sloshing scenario at the group scale. 



A.l. NGC 5044 

The dynamical status of the NGC 5044 group (z=0.009) as revealed by the pres ence of cold fro nts and sloshing features in its 
X-ray emitting gas has been discussed extensively using XMM and Chandra data ( [Gastaldello et a l. 2009 ; David et al. 2009 , and 
references therein). Here we report additional features in the X-ray morphology in this group, one of the brightest in the X-ray, 
which further support the interpretation of the sloshing scenario and the identification of t he perturber. 



In the left panel of FigfT5]we report the MOS image of the group already reported in 'Gasta ldello et al^ ( [2009 ) with a better 
color coding which highlignts the overall extension towards the East direction of the X-ray emission in the 6' -10' (66-110 kpc 
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Figure 15. Left panel: 0.5-2 keV XMM MOS image of the galaxy group NGC 5044. Point sources have been removed using the CIAO task dmf ilth. The 
excess of emission in the East direction is clearly seen. Right panel: Map of the residuals between the 0.5-2 keV MOS image of the left panel and the best fit 
beta model for the data smoothed on a 20^^ scale. The center of NGC 5044 is shown by the red cross and the position and extent of the cold fronts discussed 
in Oastaldello et al.j 2009> by the red arcs. Another relevant feature, the Eastern excess, and the sectors used for extraction of the surface brightness profiles 
discussed m the text are also shown. 
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Figure 16. Surface brightness profiles in the angular sectors shown in FigfT5] with position angle measured counterclockwise from the North direction. The 
surface brightness excess in the East can be seen in the 6^-10^ radial range. 
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Figure 17. DSS image of the central region of the NGC 5044 group. NGC 5044 and NGC 5054 are indicated. 
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Figure 18. DSS image of the galaxy NGC 5054, together with the companion N134. 



at the redshift of NGC 5044 for the cosmology adopted in this paper) radial range. In the right panel of Fig 15 we present 
the surface brightness residual map between the XMM data and th e best fit one-dimensiona l /3-model describing the surface 
brightness radial profile. The cold fronts discovered and discussed in I Gastaldello et al.] ( |2009| ) are shown and the Eastern excess 
clearly stands out in this residual map. The surface brightness profile in Fig|16|show a factor of ~ 4 enhancement compared to 
the symmetrical sector in the West. The presence of almost concentric arcs on opposite sides of the group core and the large scale 
asymmetry points towards the orbital plane of the perturber almost parallel to the LOS (e.g. Roediger et al. 2012b, and references 
therein). The SE (small scale) spiral feature extending from the SE cold front discussed in , David et al., ( ,2009, ) seems a brightness 
enhancement which connects the inner surface brightness excess marked by the cold f ront to the large scale SE excess. The same 
argument can be made for the NW small spiral, also discussed in |David et al.| ( [2009 ^, which seems to connect the core with the 
NW concentric arc marked by the outer cold front at 67 kpc. These small scale features are qualitatively seen in simulations and 
can further constrain the geometry of the merging. 

NGC 5044 is also from an optical point of vi ew an ou tstanding object in the local Universe for its notable po pulation of dwarf 
galaxies (eg. [Ferguson & San dage 1990; Mat hews et al 2004; Fa ltenbacher & Mathews |2005t|Cellone & Buzzo ni 2005; |Mendel| 



et al.||2008[|2009t|Buzzoni et al.||2012i ). |Mendel et al.[ (2008) using new, deep spectroscopic observations in conjunction with 
archival velocity information provided a sample of 111 spectroscopically confirmed memb ers. We used that sample to gain 
further information on the dynamical status of the NGC 5044 group. As already reported by "Mende l et al.| ( [2QQ8| ) the bi weight 
location and scale estimators give a mean group velocity of 2577 ±33 km s~^ and a velocity dispersion of 33 1 =b 26 km s~^ . Using 
a velocity for the NGC 5044 galaxy of 2733 ±51 km s~^ we calculate that its peculiar velocity is 156 ± 61 (just adding errors in 
quadrature) corresponds to a Z-score of 0.473 ±0.098 therefore the velocity offset is significant at the 4.8cr level. 

As for the candidate perturber, Me ndel et al.| ( |2008| ) reported a low-mass sub structure detected by means of the A test with 6 
members '^1.4 Mpc from the group center to the NE. However as suggested by David et al. (2009) this subgroup is too far away 
in the outskirts of the group to be a likely candidate, also because it has a very small component of the velocity in t he LOS (the 



biweight mean velocity of the subgroup is 2819 kms~^). A more likely candidate suggested by |David et al.| ( |2009| ) is the giant 
perturbed spiral NGC 5054 (see Fig 18 ): it is one of the brightest m embers of the group and it has a stellar mass of 2.5 x 10^^ M© 
comprising 11% of the total stellar mass of th e group dBuzzoni e t al. 2012). It is at a projected distance of 300 kpc from NGC 
5054 and it has a velocity of 1658 kms~^ ( [Mendel et al.|2008 ) corresponding to a peculiar velocity of -919 ± 61 kms~^ with 
respect to the group velocity. At a projected distance of 2.6 ' (29 kpc) from NGC 5054 lies the Sm/Im galaxy N134, one of the 
bluest dwarf galaxies of the group, partially overlapping with the north arm; their interaction is likely the cause of the intense star 
formation activity in N134 and whether or not the interaction is also responsible for the pecu liar arm m orphology of NGC 5054 
is a matter of debate ( [Buzzoni et al.|2012| and references therein). However we reinforce the David et al. (2009) suggestion that 
NGC 5054 is the likely perturber given also the X-ray evidence presented in Fig[T5]and suggest that the peculiar morphology of 
NGC 5054 is mainly caused by interaction with the global tidal field of NGC 5044 group halo. 

A.2. NGC 5846 

NGC 5846 is a nearby group (z=0.0056) wh ere evidence fo r gas sloshing with multiple cold fronts has been presented using 
a deep (120 ks) Chandra exposure in Machac ek et al. ( 2011 ). Numerous ring-like structures are present, in particular the SW 
feature was already pointed out in the XMM observation presented in |Finoguenov et al.| ( [2006|| ). In the left panel of Fig [19] we 
report the MOS image of the group (taken from the XMM observation obsID 0021540501) and in the right panel the surface 
brightness residual map between the XMM data and the best fit on e-dimensional /3-model describing the surface brightness radial 
profile. The out er cold fronts discussed in Machacek et al.[ ( [201 1\ are shown: The residual map is in agreement with the Chandra 
one reported in [Machacek et al.[ ( [2011| see their Fig.4) with a less pronounced tail like feature given the shorter (15 ks) XMM 
e xposure. 

Mahdavi et al. ( [2005 ) performed a photometric and spectroscopic survey of a 10 deg^ surrounding NGC 5846 and together 



with ROSAT, XMM, and ASCA data they suggested that NGC 5846 together with the other giant elliptical NGC 5813 are the 
dominant components of two distinct subgroups separated by 600 kpc in projection and embedded in a larger 1.6 Mpc diameter 
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Figure 19. Left panel: 0.5-2 keV XMM MOS image of the galaxy group NGC 5846. Point sources have been removed using the CIAO task dmf ilth. The 
excess of emission in the East direction is clearly seen Right panel: Map of the residuals between the 0.5-2 keV MOS image of the left panel and the best fit beta 
model for the data smoothed on a 20^^ scale. The center of NGC 5846 is shown by the red cross and the position and extent of the outer cold fronts discussed in 
[Machacek et al._ ( 201 1 J by the red arcs. 
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Figure 20. DSS image of the central region of the NGC 5846 group. NGC 5846 and NGC 5850 are indicated. 

dark matt er halo. Both the N GC 5813 and NGC 5846 galaxy groups are in the ROSAT-ESO Flux-Limited X-Ray (REFLEX) 
catalog (B ohringer et al.|2004| ). Giv en the presence of the su bgroup associated with NGC 5813 we restricted the spectroscopically 
confirmed sample of 83 galaxies of |Mahdavi et al.| ( [2005| ) to the members within 400 kpc (0.9 deg) of NGC 5846 resulting in 
a sample of 39 objects. The bi weight location and scale estimators give a mean subgroup velocity of 1833 ±7 8 kms~^ and 
a velocity dispersion of 386 ±41 kms~^ Using a velocity for the NGC 5846 galaxy of 1714 ±5 kms~^ (.Trager et al.lbOOO l 
we calculate that its peculiar velocity is -1 19 ± 78 corresponding to a Z-score of -0.339 ± 0.206 therefore the velocity offset is 

marginally significant a t the 1.6a lev el. 

As suggested by Mac hacek et al.| ( [2011| ) a likely perturber is the big spiral galaxy NGC 58 50 (see Fig|20l) which is loca ted 
10' (71 kpc) in projection to the East of NGC 5846 and it has a velocity of 2556 ±4 kms"^ ( de Vaucouleu rs et al.||l99lt . It 
therefore has a relative velocity of 842 ± 6 km s~^ with respect to NGC 5846 and of 723 ±78 km s~^ with respect to the mean 
group velocity. The optical image reveal numerous morphological peculiarities, particularly in the spiral arms and structural 
asymmetries in the HI gas in the direction of NGC 5846 are even more pron ounced (iHig don et al. 1998 ). These features already 
suggested a recent (< 200 Myr) tidal interaction between the two galaxies ( JHigdon eTal. 1998) : the new element suggested by 
the X-ray evidence is that the interaction has been between NGC 5850 and the overall mass of the group halo interior to the orbit 
of NGC 5850. 

B. REDSHIFT UNCERTAINTIES IN THE CATALOGUE OF MEMBERS OF THE IC 1860 GROUP 

Following lOwers et al. (2011b)) we performed an independent check of the measured velocities and their uncertainties for 
the galaxies in the IC 1860 group by cross-correlating the catalogues from different sources. We used as referen ce the more 
homogeneous 2d F sa mple which consists o f 74 objects a nd correlated it with the 6d FGS, the , Stein, ( ,1996, ), the [Dressier & 
Shectman ( 1988 ), the da Costa et al.| ( p^98] ) and the RC3 ( de Vaucouleurs et aLJI 1 99 1 1 ) redshifts catalogues. The results of the 
comparisons are shown in Table |4| 

As can be seen there are no significant offsets from a zero mean for all the redshift differences between the various catalogues. 
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Table 4 

Comparison of 2dF redshift measurements with literature values. 



Literature source 


'^match 


Acz 


cr(Acz) 


Literature redshift 


Quadrature difference 


Ratio 






kms-i 




uncertainty 


kms-i 


kms-i 




6dF 


10 
9 

23 


-17 ±20 
-7 ±18 
-61 ±22 


52±11 
60 ±13 
90 ±13 


31 

23 
45 




42 ±14 

55 ±14 
78 ±15 


0.35 


lSteinVl996^ 


0.26 


Dressier & ^nectman J1988| 


0.50 


11 


-30 ±29 


75 ±22 


43 




61 ±27 


0.48 


Ida Costa etal.|(|1998) 


14 


-46 ±19 


56 ±15 


37 




42 ±18 


0.41 



Note. — Columns are as follows: (1) Origin of external redshift catalog for comparison; (2) number of redshifts in common with the 2dF, Nmatch\ 0) the 
mean redshift difference, Acz = czidF - czext. , as determined by the biweight estimator; (4) dispersion in the Acz distribution, cr(Acz), determined with the 
biweight estimator; (5) the redshift uncertainty for the external redshift measurements; (6) the quadrature difference between a(Acz) and the external redshift 
error which gives an external estimate of the uncertainty on the 2dF measurements; (7) the ratio of the external uncertainty measurement to the median uncertainty 
measurement of the 2dF sample which is 89 km s~^ . The errors for Acz and a(Acz) are obtained as one standard deviation of the distribution obtained by 10,000 
bootstrap resampHngs of the data. 

ensuring that also the few external redshifts are consistent with the 2dF ones. 

C. THE NARROW ANGLE TAIL RADIO GALAXY IC 1858 





Figure 21. Left panel: GMRT image (contours and gray scale) at 325 MHz of the NAT associated with IC 1858. The restoring beam is 44^^ x 40^^ in p. a. 0°. 
The rms noise level is 1 mJy beam~^. Contours are spaced by a factor of 2 starting from 3 mJy beam~^. The inset shows the GM/?r full-resolution (14^^ x 9", 
-22°) contours overlaid on the DSS image of IC 1858. Contours are spaced by a factor of 2 starting from 1.2 mJy beam~^. The rms noise level is 0.36 mJy 
beam~^ . Right panel: 1.4 NVSS image of IC 1858. The restoring beam is 45^^ x 45^^ in p.a. 0° . The rms noise level is 0.5 mJy beam~^ . Contours are spaced by 
a factor of 2 starting from 1 mJy beam~^ . 



Radio galaxies located in dense environments often show complex and prominently distorted radio structures. A common 
morphology is represented by tailed radio galaxies, i.e., double sources whose jets and lobes are bent. Sources whose jets and 



lobes form a small angle are referred to as narrow-angle-tail radio galaxies (NAT, or U-shaped; IFeretti & Venturi|[2QQ2| and 
references therein). Many NAT galaxies resi de in clusters, however an increasing number are found in lower mass environments 
like group of galaxies (Free land et al. 12008) and references therein). Here we report on the GMRT data showing a typical NAT 
morphology for the IC 1858 galaxy (see left panel of Fig 21 ). The total flux at 325 MHz is 880 ± 70 mJy and it has a largest 
angular size of '^ 11' corresponding at ~ 300 kpc at the reashift of the source. It is also detected in the NVSS with a flux of 
~ 180 mJy in the same region of the 325 MHz emission (see right panel of Fig|2T]) implying a spectral index of a ~ 1.1 typical 
for such extended sources. Its radio power at 1.4 GHz is 2.0 x 10^^ W Hz~\ again typical for NAT sources. It is at a projected 
distance of 216 kpc from IC 1860 and it has a velocity difference of -697 ± 99 km s~^ with respect to the mean group velocity. 
Hence there is a substantial component of the velocity of the galaxy along the LOS; the overall bending of the radio jets points 
for a component in the plane of the sky pointing in the N-NE direction. The galaxy could therefore just falling in for the first time 
but also, given the it probably has already time to experience the action of the group environment transforming its morphology 
into SO, having already experience its first core passage. 



